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Abstract 

We consider a model with a small explicit breaking of a global symmetry, as suggested by 
gravitational arguments. Our model has one scalar field transforming under a non-anomalous 
U{\) symmetry, and coupled to matter and to gauge bosons. The spontaneous breaking of the 
explicitly broken symmetry gives rise to a massive pseudo-Goldstone boson. We analyze thermal 
and non-thermal production of this particle in the early universe, and perform a systematic study 
of astrophysical and cosmological constraints on its properties. We find that for very suppressed 
explicit breaking the pseudo-Goldstone boson is a cold dark matter candidate. 

PACS numbers: ll.30.Qc, 14.80.Mz, 95.35.+d 
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I. INTRODUCTION 



It is generally believed that there is new physics beyond the standard model of particle 
physics. At higher energies, new structures should become observable. Among them, there 
will probably be new global symmetries that are not manifest at low energies. It is usually 
assumed that symmetries would be restored at the high temperatures and densities of the 
early universe. 

However, the restoration of global symmetries might be not completely exact, since 
Planck-scale physics is believed to break them explicitly. This feature comes from the fact 
that black holes do not have defined global charges and, consequently, in a scattering process 
with black holes, global charges of the symmetry would not be conserved Q]. Wormholes 
provide explicit mechanisms of such non- conservation P|. 

In the present article, we are concerned with the case that the high temperature phase is 
only approximately symmetric. The breaking will be explicit, albeit small. In the process 
of spontaneous symmetry breaking (SSB), pseudo-Goldstone bosons (PGBs) with a small 
mass appear. We explore the cosmological consequences of such particle species, in a simple 
model that exhibits the main physical features we would like to study. 

The model has a (complex) scalar field transforming under a global, non- anomalous, 
f/(l) symmetry. We do not need to specify which quantum number generates the symmetry; 
it might be B-L, or a family U{1) symmetry, etc. We assume that the potential energy for 
\1/ has a symmetric term and a symmetry-breaking term 



where A is a coupling and v is the energy scale of the SSB. This part of the potential, as well 
as the kinetic term are invariant under the U{1) global transformation ^ — > e*°\l/. 

Without any clue about the precise mechanism that generates Vnon-sym, we work in an 
effective theory framework, where operators of order higher than four break explicitly the 
global symmetry. The operators would be generated at the Planck scale Mp = 1.2 x 10^^ 
GeV and are to be used at energies below Mp. They are multiplied by inverse powers of 
Mp, so that when Mp — » oo the new effects vanish. 




non—sym 



(1) 



The symmetric part of the potential is 



(2) 
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The U{1) global symmetry is preserved by vl/*vl> = |v]/|2 jg violated by a single factor 
\E'. So, the simplest new operator will contain a factor 

Vnon-sym = 'dj^l'^r (^e"^^ + ^"'c'') (3) 

with an integer n > 4. The coupling in ^ is in principle complex, so that we write it 
as g e~^^ with g real. We will consider that Vnon-sym is small enough so that it may be 
considered as a perturbation of Kj/m- Even if Q is already suppressed by powers of the 
small factor v/Mp, we will assume g small. In fact, after our phenomenological study we 
will see that g must be tiny. 

To study the modifications that the small explicit symmetry breaking term induces in 
the SSB process, we use 

^ = (p + v)e'^^'' (4) 
with new real fields p{x) and the PGB 9{x). Introducing (jH) in Q we get 

Vnon-sym = -2 g {w~^ ^'^^ ~ ~^ 

The dots refer to terms where p{x) is present. We see from (jS)) that there is a unique vacuum 
state, with < 6 >= 6 v. To simplify, we redefine 6' = 6 — 6v, and drop the prime, so that 
the minimum is now at < >= 0. From (0) we easily obtain the 6 particle mass 

(\ n-l 

Although for the sake of generality we keep the n-dependence in ©, when discussing 
numerics and in the figures we particularize to the simplest case n = 4, with the operator 
in (jHl) of dimension five. We will discuss in Section |3 what happens for n > 4. 

To fully specify our model, we finally write the couplings of the PGB Q to other particles. 
We have the usual derivative couplings to fermions j^|, 

I^^Sf = ^ (9,6) /V75/ = 9eff- h,f (7) 

We have no reason to make g' very different from 0{1). To have less parameters we set 
g' = 1 for all fermions and discuss in a final section about this assumption. Then 

9eff = ^ (8) 
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The PGB 6 couples to two photons through a loop, 

^d-y-y = o 5'6»77 ^'^"''^ FfiuFpr 9 (9) 
O 

where the effective coupling is ge^^ = In the same way, there is a coupling to two gluons, 
'^999 ~ I dsgg ^^'^^^^'fiu^pr ^ with gogg = Couphngs of 9 to the weak gauge bosons do 
not play a relevant role in our study. 

We have now the model defined. It has two free parameters: g, which is the strength 
of the explicit symmetry breaking term in (jH)), and v, which is the energy scale of the SSB 
appearing in (j21). 

At this point we would like to comment on the similarities and differences between our 9- 

n 

particle and the axion There are similitudes because of their alike origin. For example the 
form of the coupling ((Tj) to matter and © to photons is entirely analogous. A consequence is 
that we can borrow the supernova constraints on axions to constrain our PGB. Another 
analogy is that axions are produced in the early universe by the misalignment mechanism 
0] and by string decay [7] and 6'-particles can also arise in both ways. 

However, the fact that the 6'-particle gets its mass just after the SSB while axions become 
massive at the QCD scale introduces important differences. While in the axion case clearly 
only the angular oscillations matter and always one ends up with axion creation, in our 
model we need to work in detail the coupled evolution equations of the radial and angular 
part. We will need to establish when there are and when there are not angular oscillations 
leading to PGB creation. 

Also, from the phenomenological point of view, the axion model is a one-parameter model 
while our model has two parameters v and g. As we will see, this extra freedom makes the 
supernova constraints to allow a relatively massive ^-particle, and we will need to investigate 
the consequences of the decay in the early universe. There is no analogous study for axions, 
simply because the invisible axion is stable, in practical terms. 

Let us also mention other previous work on PGBs. In ^8], the authors consider the SSB 
of B-L with explicit gravitational breaking so that they obtain a massive majoron. An 
important difference between their work and ours is that we consider non-thermal produc- 
tion, which is crucial for our conclusions about the PGB being a dark matter candidate. 
Also, there is previous work on explicit breaking of global symmetries 0|, and specifically on 
Planck-scale breaking [l^. Cosmological consequences of some classes of PGBs are discussed 



in Finally, let us mention a recent paper where a massive majoron is considered in 
;lie context of a supersymmetric singlet majoron model An objective of the authors of 
1^ is to get hybrid inflation. Compared to our work, they consider other phenomenological 
consequences that the ones we study. 

The article is organized as follows. We first discuss the particle properties of 9. In 
Section IIIH we work out the cosmological evolution of the fields, focusing in thermal and 
non-thermal production, and discuss the relic density of PGBs. The astrophysical and 
cosmological bounds, and the consequences of the decay, are presented in Section IIVI A 
final section discusses the conclusions of our work. 

II. e - MASS AND LIFETIME 

We have deduced the expression (jH)) that gives the mass of the ^-particle as a function 
of V and g. In terms of these two parameters, we plot in Fig^ lines of constant mass, for 
masses corresponding to the thresholds of electron, muon, proton, bottom, and top final 
decay, = 2me, 2m^, 2mp, 2mf,, and 2mt. 

The lifetime r of the ^-particle depends on its mass and on the effective couplings. A 
channel that is always open is ^ — 77. The corresponding width is 

r(»^„) = ^ = 4,^ (10) 

where B^^ is the branching ratio. When me > 2m f, the decay 6 —>■ ff is allowed and has a 
width 



2 



where (3 = ^ 1 - 4 (j^ 

For masses me < 2me the only available decay is into photons. When we move to higher 
6 masses, as soon as me > 2me, we have the decay 6 e^e^ . Actually, then the 9 e^e~ 
channel dominates the decay because 6' ^ 77 goes through one loop. However, r(^ — > 77) oc 
whereas r(^ ff) oc m^. By increasing me we reach the value me — 150 me, where 
r(^^ ^ 77) ^ r(^^ e+e~). For higher mg the channel 6 ^ '-fj dominates again. If we 
continue increasing the 6 mass and cross the threshold 2m^, the decay 6 —>■ fj^^fi~ dominates 
over 6 e^e~ and 9 77, because Qefj is proportional to the fermion mass mf. This 
would be true until me — 150 m^, but before, the channel 9 ^ pp opens, and so on. For 



larger masses, each time a threshold 2m j opens up, 6* — > // happens to be the dominant 
decay mode. 

Now we can identify in FigI21 the regions with different lifetimes. We start with the 
stability region r > — 4 x 10^'^ s the universe lifetime, a crucial region for the dark 
matter issue. When < 2me, r = to along the dashed line in Fig^ until f ~ 4 x 10^^ 
GeV and along the line rng = 2me for higher v. Below these lines, relic 6 particles would 
have survived until now; in practical terms, for values of g and v in that region, we consider 
6' as a stable particle. Above the lines, r < tg and the particle is unstable. For future 
use, we define the ranges: (l)to > t > lO^^s; (2) 10^^ s > r > lO^s; (3) lO^s > r > lOS; 
(4) lOS > r > 10^ s; (5) lOS > r > 300 s; and (6) 300 s > r > 1 s. 

III. COSMOLOGICAL PRODUCTION AND PGB DENSITY 

^-particles can be produced in the early universe by different mechanisms. There are 
non-thermal ones, like production associated with the 6 field oscillations, and production 
coming from the decay of cosmic strings produced in the SSB. Also, there could be thermal 
production. In this section we consider these production mechanisms and the PGB density 
resulting from them. 

Let us begin with the 6 field oscillation production. The expansion of the universe is 
characterized by the Hubble expansion rate H, 



These relations between H and the time and temperature of the universe, t and T, are 
valid in the radiation era. For the period of interest, we take the relativistic degrees of 
freedom g^^ = 106.75 In the evolution of the universe, phase transitions occur when 

a symmetry is broken. When the universe cools down from high temperatures there is a 
moment when the potential starts changing its shape, and will have displaced minima. This 
happens around a critical temperature ~ v and a corresponding critical time tcr- 

In our model, the explicit symmetry breaking is very small and the evolution equations 
of p and 9 are approximately decoupled. The temporal development leading p and 9 to the 
minimum occurs in different time scales, since the gradient in the radial direction is much 
greater than that in the angular one. The evolution happens in a two-step process. First 




(12) 
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p goes very quickly towards the value < p >= 0, and oscillates around it with decreasing 
amplitude (particle decay and the expansion of the universe work as a friction). Shortly 
after, p will be practically at its minimum. In this first step where p evolves, 9 does not 
practically change, maintaining its initial value Oq, which in general is not the minimum, 
i.e., ^0 7^ 0. 9 evolves towards the minimum once the first step is completely over. This is 
the second step where the field 9 oscillates around < ^ >= 0. 

We have confirmed these claims with a complete numerical treatment of the full evolution 
equations, using the effective potential taking into account finite temperature effects. We 
summarize our results in Appendix 1X1 

Hence, we can write the following p-independent evolution equation for 9 

9{t) + 3H9{t) +771^9 = (13) 

Here overdot means d/dt. We do not introduce a decay term in the motion equation 
and will treat the 9 decay in the usual way. We can do it because 9 has a lifetime r that is 
much greater than the period of oscillations. 

Oscillations of the 9 field start at a temperature Tosc such that 3H{Tosc) ~ f^e and this 
will be much later than the period of p evolution provided Tosc ^ ^cr ~ v. This is equivalent 
to 

g^/^ < 3 X 10^^ f j (14) 

We shall see that the phenomenologically viable and interesting values of g are very tiny, so 
that we can assume that the inequality (jl4j) is fulfilled. 

Since the 9 oscillations are decoupled from p oscillations, 9 production is equivalent to 
the misalignment production for the axion and we can use those results The equation 
of motion p3|) leads to coherent field oscillations that correspond to non-relativistic matter 
and the coherent field energy corresponds to a condensate of non-relativistic 9 particles. 

We define posc and Uosc as the energy density, and number density respectively, of the 
PGBs coming from the 9 oscillations. The energy stored initially is 

Pose = rne Uosc - ^ rnl9o (15) 

The initial angle is unknown, 9o/v G [0, 27r), so that we expect it to be of order one, 9q/v ~ 1. 
In the following we will set 9q = v. Barring an unnatural fine tuned extremely close to 0, 
other choices of 9q would lead essentially to the same conclusions we reach. 



Let us now consider production of PGBs by cosmic string decays. When our U{1) symme- 
try is spontaneously broken, a network of (global) strings is formed [3| ■ These strings evolve 
in the expanding universe and finally decay into PGBs, at a time tstr given by niQtstT ~ 1, 
which is of the order of to^c- The same issue has been extensively studied for the axions and 



we can borrow the results 
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17l |. Unfortunately, there are different calculations that do 
not agree among them. To be conservative, we take the least restrictive result, namely, the 
one giving less particle production. That is ^| 

ristritstr) ~ — (16) 

This is of the same order of magnitude as the PGB density at time tstr due to 6'-oscillation 
production. 

Thus, the total number density of non-thermally produced PGBs is 

IT'non—th "^osc ~l~ "^str (l''^) 

Let us now consider the thermal production of PGBs. Any species that couples to the 
particles present in the early universe and has a production rate F larger than the Hubble 
expansion rate H during a certain period will be thermally produced. Whether such a period 



exists or not has been investigated for the axion 



m 



3- 



We adapt the axion results, taking 



into account the different magnitude of the gluon-gluon vertex, and conclude that for 

< 7.2 X 10^2 GeV = Vth (18) 

there is always thermal production of 9 in the early universe. However, for larger values of 
f , 9 interacts so weakly that T < H always, oi T > H only for such a brief period of time 
that in practice there is no thermal production. We denote this region by (I) in FigOl 

When V < Vth) the 9 species actually interacts with the plasma in the following range of 
temperature T 

V>T > = Tend (19) 

- - 1.8 X IQi^GeV ^ ^ 



and a thermal population of 9 is created with a number density 

nti^ = ^T^ 0.12 (20) 



We should now reconsider the fate of the non-thermal PGBs in (^7j), since a thermaliza- 
tion period fT^ may be at work. For this discussion we shall use the fact that 6'-oscillations 
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and string decay start at about the same temperature, Tosc ~ Tstr- Let denote this com- 
mon temperature where non-thermal production starts by T*. When v < Vth, we have to 
distinguish between two possible cases. First, if T^.. < T^nd, non-thermal 6 are born after 
the thermalization period is over. We end up having both a thermal rith and a non-thermal 
nnon-th densities, and a total density given by the sum ng = rith + nnon-th- This first case cor- 
responds to values of v and g indicated in Figl^lby regions (II) and (III); in (II) Unon-th > nth, 
while in (III) rith > Unon-th- The second case corresponds to having T* > Tend] if this is the 
case, non-thermal PGBs will be in contact with the thermal bath and will consequently ther- 
malize. Indeed, independently of details, we end up the period corresponding to (fTTHl with 
= nth given in (j2(J|) . The region (IV) in Fig.Q is where such complete thermalization 
happens. 

A word of caution is now in order. When T* > f , radial and angular oscillations are not 
decoupled. The analysis is not as simple as we presented before: the oscillations cannot 
be approximated by harmonic ones, and depend on initial conditions. However, if f < vth, 
the PGBs get thermalized in any case, so that we have not to worry about the non-thermal 
production details. For v > Vth (and still T^, > v), since there is no thermalization, the final 
density is ng = Unon-tk- Admittedly, we have no general formula for the density in this case, 
represented by (V) in Fig|21 In this region, 9 lifetimes are so small that the particle does 
not play any cosmological role at all. 

Summarizing, a number density of 9 particles always appears in the early universe. In 
regions (II), (HI), and (IV), the thermal density is rithiT = To) at decoupling temperature 
Td (the temperature where F = H). In regions (I), (II), and (III) there is also a non-thermal 
density nnon-th at formation temperature T = T^. 

As usual, the expansion of the universe dilutes these densities. Eventually, if the particle 
is unstable it will decay. The issue of unstable 9 and effects of the decay will be worked out 
in Section IIV| and here we focus on stable 9 particles, since certainly they would be dark 
matter. To calculate the relic density of 9 today, apart from the expansion effect, one has to 
take into account the transferred entropy to the thermal bath, due to particle-antiparticle 
annihilations. We follow the standard procedure (see ^^), and find that the number density 
today is 
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where Ti = To for rith and Ti = T^: for Unon-th, and Tq = 2.7-?^ is the photon temperature 
today. We take g^s{Ti) = 106.75 and g^s{To) = 43/11. 

With all that, we can calculate the cosmological density of 6 particles that we would have 
today, n0{to), and the corresponding energy density normalized to the critical density Pc, 

= ^^^^^ (22) 

Pc 

with = 0.5 X lO-^GeVcm-3 Q. ^ 

All cosmological data available lead to a dark matter contribution Qdai — 0.3 [14]. In 
FiglUwe plot the line Qg = 0.3. We will see in Section HVl that part of this line is excluded. 
Clearly, values of g and v along the non-excluded part of the line would imply the density 
needed to fit the cosmological observations. Our 6 particle is a dark matter candidate 
provided g and v are on the line, or not far, in such a way that is still a substantial 
fraction of 0.3. We stress that in the region of interest for dark matter it is the non-thermal 
production that dominates (see FiglHl). Thus, the produced PGBs are non-relativistic and 
would be cold dark matter. Values for Qg much smaller than ~ 0.3 are allowed but not 
cosmologically interesting. On the other side of the line, but still in the stability region, we 
have the excluded region Qg > 0.3. 

IV. ASTROPHYSICAL AND COSMOLOGICAL CONSTRAINTS 

PGBs would be emitted from the hot stellar cores since nucleons, electrons and photons 
initiate reactions where 6 is produced. Provided 6 escapes the star, the emission leads to a 
novel energy loss channel, which is constrained by stellar evolution observations using what 
is nowadays a standard argument Isf . The limits we find for v are similar to the limits found 

n 

for the Peccei-Quinn scale i2(j\ because the coupling of 6 to matter is the same of the axion. 
We get the bound 

v>3.3x 10^ GeV (23) 

that is shown in FigEl 

We see that the astrophysical constraint is on v, but not on ^f, so that we may have 
relatively massive PGBs, while in the axion case this leads to < O(10~^eV), i.e., to 
forbid the existence of relatively massive axions. This makes the PGB in our model an 
unstable particle, at variance with the invisible axion. 
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Let us now investigate the region where 9 is unstable, see Fig0J The cosmological effects 
of the decay will depend on the lifetime r, mass m^, and number density ne at the moment 
of the decay. These properties depend on the parameters g and v. Observational data will 
help us to constrain different regions of the v parameter space. We will use several pieces 
of data, depending on when the 6 decay occurs. To do this systematically, it is convenient 
to distinguish between the following time ranges, plotted in Fig|21 
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We do not consider times earlier than tsBN since there are no observational data that 
give us any constraint at that times. The physical meaning of the chosen times are 

- tree IS the time of recombination, i.e., the moment at which photons last scattered with 
matter, when free electrons present in the cosmic plasma were bounded into atoms; 

- tc is the time at which the rate of Compton scattering ej — >■ 67 becomes too slow to 
keep kinetic equilibrium between photons and electrons; 

- tth is the time at which the double- Compton scattering 67 677 decouples; 

- ^He is the time at which the energy of the Cosmic Background Radiation (CBR) is low 
enough to permit photons of energy ~ 10 MeV to photodissociate ^He nuclei; 

- tendBBN is the time at which the Big Bang Nucleosynthesis (BBN) finishes; 

- ^BBN is the time at which BBN begins. 

We can constrain the parameters v and g because the photon spectrum would be distorted 
by 6 decay products, and the abundances of the light elements would be altered. We will 
discuss this in the next two subsections. 
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Effects on the photon spectrum 



The products of the 9 decay can be either photons or fermions. When the decay products 
are fermions it is important to note that if the number density of the created fermions is 
higher than the photon density of the CBR, these fermions immediately annihilate into 
photons, since we are considering decays produced after the e^e~ annihilation in the early 
universe. The distinction between the two channels will be important only when the 6'-decay 
occurs in the time ranges 1 and 2, since in the other regions photons and fermions thermalize. 

When T is in the range 1, it is convenient to distinguish between two regions (see FigsQ 
and 121) • In the region where itlq < 2me, 9 decays only into photons. The products of the 
other region are mainly e~ and e"*", but the values of v and g in this region are such that 
the number density of the created fermions is higher than the CBR photon density. Then, 
in all the time range 1, the final products are photons, but for nig > 2 me we use r for the 
9 e^e~ decay. Photons produced in this range stream freely and contribute to the diffuse 
photon background of the universe. One can compute the present energy flux of photons 
per energy and solid angle interval coming from 6'-decay ^ 



d^J^E 1 rieito) f E 

exp 



dEdn 27T tHo \me/2 



2 1 f E 



3 T Hq \'me/2 



3/2' 



(24) 



where Hq is the Hubble constant today and n0{tQ) is the number density of 9 that we would 
have today if it had not decayed. The predicted flux is bounded by the observational 
limits [2^ on the photon background flux in the energy range of interest. This restriction 
excludes all the region in the v, g parameter space corresponding to time range 1. Taking 
into account that even for t > to a small fraction of PGBs have decayed into photons, we can 
also prohibit part of the parameter space that corresponds to stable particles with Q < 0.3 
(see Section ITTT|) . In Fig. |S]we plot this excluded region. 

In range 2, Compton scattering is not effective in maintaining kinetic equilibrium between 
the e~ and the 7 of the CBR. Depending on the precise values of g and v, 9 ^ e~^e~, 
9 or — i> 77 will be the dominant decay. The induced charged leptons, carrying 

large energy, outnumber the existing CBR electrons. Different processes now compete; one is 
scattering of these hot electrons and muons with CBR photons. Another is e^e~ and /i 
annihilations that give high energy photons, which heat CBR electrons. The last is high 
energy photons produced in the decay 9 77, that also scatter and heat CBR electrons 
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which in turn scatter with CBR photons. Whatever process is more important (it depends 
on the (7, V values), the result is a distortion of the photon spectrum, parameterized by the 
Sunyaev-Zeldovich parameter y j23|. The energy /S.E dumped to the CBR, relative to the 
energy of the CBR itself is constrained by data on CBR spectrum 

-Ay< 4.8 X 10^5 (25) 

^CBR 

where = mgUQ. The experimental bound (j25j) rules out all the g and v values that 
would imply a lifetime r in the time range 2. 

In the range 3, Compton scattering is fast enough to thermalize the products of the 9- 
decay occurring in this range. This is because even in the region where the products are 
neither photons nor e^e~, the final particles will be photons in any case. In this region 
the 67 — > 677 processes are not effective so the photon number cannot be changed. So, 
after thermalization, one obtains a Bose-Einstein CBR spectrum with a chemical potential, 
/ = [exp (E + /x)/T — 1] ^, instead of a black-body spectrum. The relation between and 



AE is 



24| 



4 C(2) C(3)^ .,_ 4 An^ 



3 C(3) C(4); ^ EcBR 3 n, ^^^^ 
CJji) is Riemann's zeta function). The parameter /i is very well constrained by CBR data 
201 that gives < 9 x 10~^. This value rules out all the v region that would give r in 
the time range 3. 

For times before the range 3, both Compton and double Compton scattering are effective, 
so the decay products thermalize with the CBR, without disturbing the black body distribu- 
tion but changing the evolution of the temperature of the thermal bath. This temperature 
variation leads to a change in the photon number, and thus to a decrease on the parameter 
rih/n^- The knowledge we have on the value of this parameter at tree and tendBBN 
allows to constrain the g^ v region that would give lifetimes r in the time ranges 4 and 
5. Although the corresponding restriction is quite poor (it is essentially AE/Ecbr ^1), yet 
it excludes the parameter space corresponding to the ranges 4 and 5. However, constraints 
from the effects of the 6'-decay on the light element abundances are much more restrictive 
in these ranges, as we will examine in the next subsection. 
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Effects on the abundances of the hght elements 



The period of primordial nucleosynthesis is the earliest epoch where we have observational 
information. Also, the theoretical predictions of the primordial yields of light elements are 
robust. The agreement with observation is considered one of the pillars of modern cosmology. 
The 9 decays, and the 9 particle itself, might modify the abundances of the light elements, 
which implies restrictions on the f , g parameters. 

First, we consider how the decays of a PGB affect the abundances of the light elements 
after they are synthesized, i.e., after tendBBN- One of the consequences of the decay is the 
production of electromagnetic showers in the radiation-dominated plasma, initiated by the 
decay products. As a result, photons of energy ~10 MeV scatter and photodissociate light 
elements. This scattering occurs after tne, because before tne, the collision of these photons 
with the GBR is more probable than with the light e ements. In time range 4, the photodis- 
sociated element is deuterium (if uiq > 10 MeV) j25| while in time ranges 3 and 2 there is 
helium photodissociation (if > 100 MeV), with the subsequent production of deuterium. 
Observational data for the abundance of deuterium constraint all these processes. When 9 
decays into quarks which hadronize subsequently, hadronic showers can also be produced 
(if mg ~ 1 GeV). These processes dissociate ^He before twe, overproducing deuterium and 
lithium. All these constraints, that are summarized in j25|, exclude all the g and v values 
that give r in the time ranges 2,3,4 and 5, provided mass conditions are fulfilled for each case. 
However, values from ranges 2,3,4 and 5 that do not satisfy the proper mass restrictions, 
are anyway excluded by the constraints we considered in the former subsection. 

The other effect on the light element abundances arises because the 9 particle would 
modify the BEN predictions. The presence of 9 and, more important, the presence of the 
relativistic products of the ^-decay, accelerate the expansion rate of the universe in the 
relevant BEN period and modify the synthesis of the light elements. The decay of the 9 
boson is also a source of entropy production, which alters the temperature evolution of the 
universe. This changes the number of photons (and the value of t]) and produces an earlier 
decoupling of neutrinos. Then the relation between and is changed, with potential 
effects on the EEN physics. It is important to note that this production of entropy never 
rises the temperature of the universe and does not lead to several successive EEN 
periods. All these effects in EEN have been implemented in the usual code, which 
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allows to constrain the quantity ^E/n^. As a result, some of the values of (7, v that would 
give r in the time range 6 are ruled out. The potential effects of hadronic showers, that we 
have mentioned earlier, also would modify the BBN results [25[, allowing us to exclude part 
of the range 6, but not all of it. 

All the results exposed in this section are summarized in FigEl There, we can see that 
after our systematic study, it turns out that all the zone of parameters that lead to lifetimes 
between about 1 s and to is ruled out. Only a small zone corresponding to range 6 is allowed. 



V. DISCUSSION AND CONCLUSIONS 



Gravitational arguments suggest that global symmetries are explicitly violated. We de- 
scribe this violation using an effective theory framework that introduces operators of order 
higher than four, suppressed by inverse powers of the Planck mass. These operators are 
considered as a perturbation to the (globally) symmetric part of the potential. 

The SSB of global symmetries with a small explicit breaking leads to PGBs: Goldstone 
bosons that have acquired a mass. Equivalent to the appearance of a mass for the boson, 
there is no longer an infinity of degenerate vacua. 

In this article we have studied the cosmology of the PGB 9 that arises in a model con- 
taining a scalar field with a potential that can be divided in a part that has a global non- 
anomalous U{1) symmetry and another part with gravitationally induced terms that are 
U{1) violating. In our analysis we let vary two parameters of the model: the SSB breaking 
scale, V, and the coupling g of the Planck- induced term in the potential. We have analyzed 
the evolution of the field towards the vacuum in a phase transition in the early universe. It 
occurs in a two-step process: first the radial part attains its minimum in a relatively short 
time and second the angular part of the field starts oscillating well after the first step is 
over. The 9 field oscillations correspond to non-relativistic matter. We have calculated the 
density of 9 particles born through this mechanism and also through the decay of cosmic 
strings created at the SSB. There might also be thermal production of PGBs in the early 
universe and also there might be thermalization of PGBs produced non-thermally; these 
issues have been fully analyzed in our article. 

A variety of arguments constrain the parameters v and g of 9. There are astrophysical 
constraints coming from energy loss arguments. There are also cosmological bounds. When 
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the particle is stable its density cannot be greater than about the critical density, otherwise 
the predicted lifetime of the universe would be too short. If the particle is unstable the 
decay products may have a cosmological impact. We have watched out for effects of the 
decay products on the CBR and on the cosmological density of light elements. 

We have considered all the above potential effects and used empirical data to put con- 
straints on g and v. We have been led to exclude the region of the g, v parameter space 
indicated in FigEl 

In Figini we see that there are two allowed regions in the plot. First, in the upper part of 
the plot there is an allowed region. It is where r < 1 s, except the tooth at values that are 
about V ~ 10^^ GeV and g ~ 10^^^ that corresponds to 1 s < r < 300 s (part of zone 6 in 
Fig|21). For a 6 that has the parameters corresponding to this first region, it will definitely 
be extremely difficult to detect the particle. Also, in any case, it will have no cosmological 
relevance. 

The second permitted zone of the figure is where 9 is dark matter, at the bottom of the 
plot. It would be an interesting cold dark matter candidate provided the values of g and 
V are not far from the solid line in FigEl There is an upper limit to the mass mg in the 
allowed region where 6' is a dark matter candidate 

me < 20 eV (27) 

A way to detect 9 would be using the experiments that try to detect axions which make 



ar coupling to two photons exists 
2^. The detection techniques use 



use of the two-photon coupling of the axion. Since a simi 
for the PGB, we would see a signal in those experiments 
coherent conversion of the axion to photons, which implies that in order that 9 would be 
detected, we should have me < 10"^ eV. 

For 9 be dark matter, we notice that the values of g have to be very small 

g < 10"^° (28) 

We do not conclude that these values are unrealistically small. Without any knowledge of 
how gravity breaks global symmetries it would be premature to argue for or against the order 
of magnitude (j2Hl)- For example, in j2^, Peccei elaborates about the explicit gravity- induced 
breaking of the Peccei- Quinn symmetry, and gives the idea that perhaps the finite size of a 
black hole when acting on microscopic processes further suppresses Planck-scale effects. 
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Apart from that, there is an easy way to get PGBs as dark matter candidates for values 
of g not as tiny as in (f^ . Notice that to work out the most simple case, we have considered 
n = 4 in Eq.Q. It suffices to consider more general potentials 

Vnon-sym = "^^^^^l^l" (^'"e"^^ + h.C. ) (29) 

with n,m integers. In the present article we made n = 4, m = 1. If we take greater values, 
we get a suppression of the symmetry breaking term due to extra factors v/Mp and we may 
allow values for g higher than the ones obtained for g. In order of magnitude, for v ~ 10^^ 
GeV, we see that taking operators of dimension n + m = 8 or 9 we have a PGB that is a 
dark matter candidate, but now with ^ ~ (9(1). In building the model we should have a 
reason for having the order of the operators starting at n + m > 5. The standard way is to 
impose additional discrete symmetries in the theory. 

Finally, we would like to comment on having put g' = 1 in We could, of course, 
maintain g' free, even we could let g' be different for each fermion, but, in our opinion, the 
introduction of extra parameters would make the physical implications of our model much 
less clear. This is why we fixed g' = 1, but now it is time to think what happens for other 
values of g'. 

The coupling g' appears accompanying a factor proportional to the mass of the fermion 
and inversely proportional to the energy breaking scale v, as expected for real and pseudo 
Goldstone bosons. When a fermion has a U{1) charge, we have no reason to expect that g' 
is much different than 0{1), but if a fermion has vanishing U{1) charge then the coupling of 
6 to this fermion may only go through loops, and consequently we have a smaller g'. In this 
case, an important change concerns the astrophysical bound. Since goNN is smaller than 
rriN/v, the bound from supernova is weakened and values of v smaller than 3.3 x 10^ GeV 
would be allowed. Another effect is that the bounds coming from the cosmological effects of 
the 6 decay are relaxed, since the lifetimes are longer when g' is smaller. However, this does 
not mean that part of the prohibited region in FigEl would be allowed. We have to take 
into account that non-thermal production is not altered when changing g' and therefore the 
bound Qe < 0.3 leads to strong restrictions in the g,v parameter space. Also, 6 thermal 
production is suppressed with smaller g' . 

We would like to show graphically what would happen for very small g', and with this 
objective we show in Fig|7|the permitted and the prohibited regions in the limit g' 0. At 
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the view of the result, we conclude that one may have PGB as a dark matter candidate for 
much larger values of g than obtained before in (|^H|). 



Acknowledgments 

We thank Mariano Quiros and Ramon Toldra for very useful discussions. We acknowl- 
edge support by the CICYT Research Project FPA2002-00648, by the EU network on Su- 
persymmetry and the Early Universe (HPRN-CT-2000-00152), and by the Departament 
d'Universitats, Recerca i Societat de la Informacio (DURSI), Project 2001SGR00188. One 
of us (G.Z.) is supported by the DURSI, under grant 2003FI 00138. 



APPENDIX A 



1. How to obtain the effective potential 

We present here in some detail how to find the effective potential that gives us a complete 
description of the physics involved in our model. Following the standard procedure js^, 
taking into account the finite temperature effects, we are led to a new contribution to Vgym, 
which is given by 

1 1 r°° r / 1 

VP = ^M^'f^"] = ^ I dxx' In [l - e-V^^'^] (Al) 

where Jb is the thermal bosonic function and [3 = 1/T, and rri^ = — |Af^ + A\1/*\E' + |AT^ 



is the effective mass. With ()A1|) . we see the behavior of the finite temperature effective 
potential. For practical applications, it is convenient to use a high temperature expansion 
of js^ written in the form 

yl^ ~ —m^T' - —m^T - -^m^ In "^^^ (A2) 
24 127r 647r2 223.63 ^ ' 

where we have neglected terms independent of the field. The effective potential must contain 
the explicit symmetry-breaking term of our model, Vnon-sym- Using for \1/ the parametriza- 
tion \1/ = 06*^/"", the expression for this term is 



kri+l 



Vnon-sym — 2(7"^^^^^ COS ( - ) (A3) 



Thus, our effective potential will be written as the sum of three terms 



Veff Vsyra ~\~ V -\- Vnon-sym (■^^) 
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2. How to find Tcr 

SSB of Veff is triggered at a critical temperature Tcr corresponding to time tcr, and 
instead of having a minimum at = 0, a local maximum appears. T^r is the temperature 
when the second derivative of Vg// at = cancels. In calculating the derivative of K//, 
we find that 

-<l>-mT) (A5) 



f{(l),T) contains all the information we need to find both T^r and the minimum of the 
potential as a function of T (writing f[(j),T) = and solving for 0). The second derivative 
of Veff at its origin is 



= /(0,T)|^=o + (0/'(0,T))|<^=o = /(0,T) (A6) 

Tcr is the solution to the equation /(O, T) = 0, and depends on A and v. In order of 
magnitude, Tcr ~ v. 



3. Field evolution 



At early times, the effective potential has a unique minimum at = 0. As the temperature 
goes down, its effects on the effective potential diminishes, and at Tcr it is spontaneously 
broken and a second order phase transition occurs. In this subsection, we present our study 
on this phase transition focusing on the separate evolutions of the radial (0) and angular 
(9) parts of and which are the implications of Vnon-sym in our effective potential. 

The evolution of the two fields mentioned above is described by these equations 

0(t) + {3H + r,)0(t) + ^ = 
9{t) + 3He{t) + %i =0 

where H = 1/ (2t) is the Hubble expansion rate of the universe, and overdot means time 
derivative. An important difference between the two equations is the appearance of in 
the first one, and this is because 9 couples derivatively to matter, see ((Zj), while couples 
as g^ffff(p- Therefore, the coupling of 9 is weak, since it is suppressed by a high-energy 
scale V, but the coupling of the radial part is not and has to be introduced in the evolution 
equation for 0. For our numerical simulations, we put F^ = m^p/SiT = V2Xv/8'n', which 
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corresponds to g^jj = 1 and consider only one species of fermions in the decay. The two 
differential equations ()A7|) are not independent of each other. They are related due to the 
fact that Veff has both and 6 dependence. However, the coupling between equations, in 
the case where g is tiny, is small and can be neglected in a first approximation. This is 
what we have done in Section lllll Here we do not neglect it since we want to do a complete 
numerical analysis and calculate the evolution when g is not small and check when the 
scenario described in Section ITTTl breaks down. It is convenient to replace the temperature 
dependence of K// with a time dependence, using relation (fT^ 



Here C = -^Mp considering g^ = 106.75 for the temperature range where we apply the 
evolution equations. To simplify calculations and graphical displays, we introduce new 
dimensionless variables 

0=-; = -; i=^; tcr = ^. (A9) 

V V tcr T^r 

With these changes, the two evolution equations are 

ki) + (I + lli:VA„)^(i) + ^/-^ = 

m + ¥w + = 

We can numerically solve the system (jA10|l for any values of interest for v and g. In FiglHl 
we show one such solution for arbitrarily chosen values for v,g, and A, for n = 4. 



4. Discussion 



In Section ITTTl we based our work on the fact that, due to the smallness of the explicit 
symmetry breaking, the field evolution towards its minimum occurs in a two-step process: 
firstly, the radial field goes quickly towards its vacuum expectation value, oscillates around 
it for a finite time till it stops due to the expansion of the universe and coupling to fermions; 
secondly, the angular field stays constant much longer and finally starts to oscillate around 
its minimum. We have proved numerically that this is so; we have checked it by solving the 
system equations ()A10|1 for a variety of values of our parameters. 

We would like to specify the upper limit on g for our model to make sense. The parameter 
g is considered to be too large when the term Vnon-sym in the effective potential starts to 
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dominate over the other ones, for ~ f . When this happens, the exphcit symmetry breaking 
is so big that, where the absolute minimum of the effective potential is supposed to be, the 
first derivative of Vejf with respect to is negative and there is no minimum at all. In 
this case, there makes no sense talking about angular oscillations around the minimum. 
Therefore, by comparing the (A/4)0^ term with Kion.-sj/m, one obtains an upper limit for g 

9<3^ (All) 

We have numerically checked that, for parameter values of interest to us, radial oscillations 
start very early and they are very rapidly damped, at a time scale much less than the time 
when ^-oscillations start. One example is plotted in FiglHl What happens is that the radial 
field oscillates for a small time around the minimum of the symmetric part of the effective 
potential and after the oscillations stop, the field stays at the minimum and evolves in time 
until temperature effects are irrelevant and the minimum stabilizes at = 1 {(f) = v). Thus, 
for values of g that satisfy ()A11|) and (fT^ . when angular oscillations start, the radial ones 
have already stopped and we must not worry about the possibility that the two oscillations 
happened at the same time. More important is to impose the condition that when radial 
oscillations start, the minimum of VJ=// be located close to the value cf) = v m order to have 
initial conditions for ^^-oscillations also of order v. With all this in mind, for v values in the 
range of interest (10'^ GeV< v < 10^^ GeV), we get to the conclusion that g must be smaller 
that about 10~^ (the number depends slightly on v and A). In particular, considering values 
of interest for to be a dark matter candidate {v ~ 10^^ GeV) and A ~ 10~^, we obtain 
an approximate limit g < 10~^. This upper limit is also roughly given by the dotted line 
represented in FiglHl which corresponds to T* ~ v. Obviously, values for g ~ 10^'^'^ and 
V ~ 10^^ GeV which we have found to be interesting to have 6' as a reasonable dark matter 
candidate of the universe, are consistent with our mechanism of producing 9 particles. 
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FIG. 1: From bottom to the top, the five sohd hnes are the hnes of constant mass, in the g,v plane, 
for 1710 = 2nie, 2m^, 2mp, 2mi,, and 2mt. The dashed line corresponds to the 9 lifetime equal to 
the universe lifetime, t = t^, when ^ 77 is the only available mode. 
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FIG. 2: Zones corresponding to some ranges for the 9 lifetime, r. At the bottom, we have the 
region of the parameter space leading to a stable 9 (r > Iq), and at the top, the zone 7, where r < 1 
s. The other regions correspond to : (1) to > t > 10^^ s; (2) 10^^ s > r > 10^ s; (3) 10^ s > r > 10^ s; 
(4) lOS > r > lOS; (5) 10"^ s > r > 300s and (6)300s > r > Is. 
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FIG. 3: We show the various regions according to the 9 production mechanism in the early 
universe. With this objective, we display the lines: T^, = v (dotted line), T^, = Tend (solid line), 
n-th = nnon-th (dashed line) and v = Vth (dot-dashed line). In (I) and (V) there is only non-thermal 
production. In (II) and (III) there are non-thermal and thermal processes that generate 0, in (II) 
non-thermal dominates while in (III), it is thermal production that dominates. In (IV) there is 
thermally and non-thermally produced 6s, but the last are thermalized. In (V), we have not a 
general formula for predicting Unon-th- 
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FIG. 4: Properties of 9 before considering the astrophysical constraints and the effects of the 
decay. Above the dashed hne of the plot we have the region of v where 9 is unstable. The solid 
line corresponds to = 0.3. The interesting dark matter region is labeled in the plot. The region 
between the solid and dashed lines is excluded because 9 is stable and Q.e > 0.3. 
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FIG. 5: The astrophysical constraints forbid the shadow region labeled SN1987A. The other 
shadow region corresponds to stable particles (r > to)) yet the contribution to the photon back- 
ground is too high and thus this region is forbidden. To help the reader, we also plot the line 
Clg = 0.3 (solid) and the border between stable and unstable PGBs (dashed). 
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FIG. 6: The prohibited region, using ah constraints we have studied, is in shadow. In white, we 
show the aUowed region. The sohd hne corresponds to ^1$ = 0.3. 
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FIG. 7: Permitted and prohibited regions in the hmit g' — > 0. The sohd Une is = 0.3 and the 
dotted hne is T^, = v. In the upper part of the dotted line we have no reUable way to calculate VLq. 
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FIG. 8: Evolution of the normalized fields 9 and as a function of r/ = \og{t/tcr) with tcr defined 
in ()A9|) . We see how (p evolves first while remains constant (a), and how 6 finally oscillates (b). 
Notice the (logarithmic) time scales. For this numerical simulation we chose: v = 10l^ A = 10"^ 
and g = 10~®. 
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